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I discuss the results of a new calculation of the charm and bottom quark masses in the quenched approximation 
and in the continuum limit of lattice QCD. The work has been done by the APE group at the "Tor Vergata" 
University making use of the step scaling method, previously introduced to deal with two scale problems, 
that allows to take the continuum limit of the lattice data. We have computed the RGI quark masses and then 
we have connected the results to the MS scheme. The continuum numbers are mf^^ = 6.73(16) GeV for the 
&-quark and m^^' — 1.681(36) GeV for the c-quark, corresponding respectively to [m^^) = 4.33(10) GeV 
and m^'^ (m^^^) = 1.319(28) GeV. The latter result, in agreement with current estimates, is for us a check of the 
method. 



1. Introduction 

Quark masses are fundamental parameters of 
the QCD Lagrangian. Their accurate knowledge 
is required in order to give quantitative predic- 
tions of fundamental processes. A direct ex- 
perimental measurement of quark masses is not 
possible because of confinement, and their de- 
termination can only be inferred from a theo- 
retical understanding of the hadron phenomenol- 
ogy. The calculations can be numerically per- 
formed through different strategies, depending 
upon the quark flavor and the available computa- 
tional facilities. Accurate non-perturbative mea- 
surements of the u, d, s and c quark masses have 
been obtained from a straightforward comparison 
of hadron spectroscopy and lattice QCD predic- 
tions. 

The situation is different for the b quark, with 
present computers capabilities. In principle, the 
b quark mass could be extracted, similarly to 
lighter flavors, by looking at the heavy-light and 
heavy-heavy meson spectrum. However, heavy- 
light mesons are characterized by the presence 
of two different scales, i.e. Aqcd, that sets the 
wavelengths of the light quark, and the heavy b- 
quark mass. Managing these two scales in a naive 
way would require a very large lattice (0(100"*) 
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points). 

The task of determining the 6-quark mass has 
been faced in literature by resorting to some ap- 
proximations of the full theory, e.g. HQET on the 
lattice 0, lattice NRQCD 0, or QCD sum rules 
j4l5j . A novel approach recently introduced is 
based on the non-perturbative renormalization of 
the static theory and its matching to QCD (jfTj . 

An alternative approach to the bottom quark 
physics, based on finite size scaling, has been 
proposed in a previous paper jHj, where it has 
been applied to the heavy-light meson decay con- 
stants (see also PTU] ). The main advantages of 
this step scaling method (SSM) are that the en- 
tire computation is performed with the relativis- 
tic QCD Lagrangian and that the continuum limit 
can be taken, avoiding the unfeasible direct cal- 
culation. In order to implement the SSM, a fi- 
nite size scheme is required, and we adopt the 
Schrodinger Functional (SF) as the most useful 
framework. 

2. Step scaling functions and HQET 

The SSM has been designed in order to deal 
with two scale problems in lattice QCD 8 and 
it has been discussed in detail in j9llj . The main 
assumption of the method is that the finite size 
effects affecting the meson masses, M^py^ have 
a mild dependence upon variations of the heavy 
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quark mass and are controllable from a numerical 
point of view. Finite size effects can be obtained 
from 
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for the pseudoscalar and vector meson masses re- 
spectively. 

To validate the hypothesis of low sensitivity 
upon the high energy scale, we can make use of 
the HQET predictions on the heavy-light meson 
masses. In the infinite volume the pseudoscalar 
and vector meson masses have the following ex- 
pansion in terms of the heavy-quark mass: 
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where X € {P,V}. Assuming the contribution 
of the corrections to be negligible, at finite 

volume one has 

Mx (m, ,mi,L)^ m^+Ax [mi , (3) 

mh 

where Kx{mi, L) depends upon the spin of the 
meson state because of the contamination of the 
excited states to the finite volume correlations. 
Using eqs. (QJ and Q we obtain the HQET 
predictions for the step scaling functions of the 
heavy-light meson masses 
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This result requires some considerations. First 
we want to stress that in the infinite heavy-quark 
mass limit the step scaling functions have to be 
exactly equal to one, axiL,mi,mh 00) = 1^. 
This represents a strong constraint for the fits 
of the heavy-quark mass dependence of the step 
scaling functions. 

The second observation concerns the number of 
terms to be considered in eq. At order 

0{l/mh) one has 
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corresponding to the static approximation in 
eq. By increasing the physical volume L, the 
difference between Ax {mi, 2L) and Ax {mi, L) de- 
creases because the two quantities have to be 
equal in the infinite volume limit, making the 
heavy-quark mass expansion of the finite volume 
effects rapidly convergent. The same arguments 
apply to the coefficient cr^-' {mi, L) that has to be 
considered when in the expansion of the meson 
masses, eq. (PJ, the order 0{l/mh) is taken into 
account. In our calculation we have performed 
the fits of the step scaling functions considering 
the 0(1/to^) term, cr^'', beyond the so called 
static evolution {SE) that retains cr^^' only. We 
will also report for comparison the fits for the SE 
that are anyway compatible. 
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Figure 1. The figure shows the continuum ex- 
trapolated step scaling functions ap{Lo) as func- 
tions of 1/ mf '^^ with two possible fits. The heavy 
extrapolations are shown only for the heavy- 
strange (Hs) set of data. 



3. Results and Discussion 

The physical numbers for the heavy meson 
spectrum have been obtained combining the re- 
sults of a small volume {Lq = 0.4 fm) calculation, 
where the simulations have been performed at the 
physical values of the heavy quark masses, with 
the results of the step scaling functions according 
to the identity: 



Mx{Loo) = Mx{Lo) <7x{Lo) (Tx(2Lo) 



(6) 



The step scaling functions at the values of the 
heavy quark masses simulated on the small vol- 
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ume have been obtained by interpolation in the 
heavy-Hght case using eq. Q while, in the heavy- 
heavy case, the results are linearly extrapolated. 
The plots of (Tp(Lo) and crp(2Lo) sls functions of 
the inverse quark mass are shown in figs. ^ and |2] 
respectively. 
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Figure 2. The figure shows the continuum ex- 
trapolated step scaling functions ap{2Lo) as func- 
tions of with two possible fits. The heavy 
extrapolations are shown only for the heavy- 
strange (Hs) set of data. 

Comparing the numerical heavy meson spec- 
trum with the experimental results |14| we have 
obtained different determinations of the &-quark 
mass, depending upon the physical state used as 
experimental input. 

Within the quenched approximation, the deter- 
minations of the quark masses coming from the 
heavy-heavy or from the heavy-light spectrum 
in principle differ because the theory does not ac- 
count for the fermion loops. The numbers we 
quote as final results are obtained by averaging 
the results for the heavy-heavy and the heavy- 
strange vector and pseudoscalar mesons and by 
keeping the typical error of a single case: 

77if"^^ = 6.73(16) GeV 
mf^imf^) = 4.33(10) GeV (7) 
for the &-quark and 

m^^^ = 1.681(36) GeV 
TOf^(mf^) = 1.319(28) GeV (8) 

for the charm. The latter results compare favor- 
ably with the results of the direct computations 

una. 



Our error estimate includes both the statistical 
error from the Monte Carlo simulation as well as 
the systematic error coming from the uncertainty 
on the lattice spacing corresponding at a given /3 
value (we have used ro = 0.5 fm Illl2ll3 j) and 
from the uncertainty on the renormalization con- 
stants. The final errors on the continuum quanti- 
ties, of the order of 2% percent for the renormal- 
ization constants and of about 1% percent for the 
scale, are added in quadrature and then linearly 
added to the statistical errors. The evolution to 
the MS scheme has been done using four-loop 
renormalization group equations. 
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